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We consider the dynamics of a 3-brane embedded in an extra-dimensional Tolman-Bondi Universe 
where the origin of space plays a special role. The embedding is chosen such that the induced matter 
distribution on the brane respects the spherical symmetry of matter in the extra dimensional space. 
The mirage cosmology on the probe brane is studied, resulting in an inhomogeneous and anisotropic 
four dimensional cosmology where the origin of space is also special. We then focus on the spatial 
geometry around the origin and show that the induced geometry, which is initially inhomogeneous 
and anisotropic, converges to an isotropic and homogeneous Friedmann-Lemaitre 4d space-time. For 

£X instance, when a 3-brane is embedded in a 5d matter dominated model, the 4d dynamics around 

the origin converge to a Friedmann-Lemaitre Universe in a radiation dominated epoch. We analyse 

f*"*) ■ this isotropisation process and show that it is a late time attractor. 
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I. INTRODUCTION 



The possible existence of extra dimensions is an exciting prospect. Ever since the first unification attempt by Kaluza 

O ' |l| and Klein [2| , and then the advent of string theory, the role of these extra dimensions in shaping four dimensional 

D ■ physics has been emphasized [3|t51 • It could range from the presence of towers of excitations in the particle spectrum 

rG \ to the appearance of scalar fields in the form of compactification moduli. In the last fifteen years and in particular 

after the discovery of fundamental branes and the building of phenomenological brane-world models [6[ such as the 

£\l ' one constructed by Randall and Sundrum [7|, |8[ , the emerging possibility of extra dimensional cosmology has been 

^- , considered thoroughly. In brane cosmology, a brane configuration comprising one or several branes are embedded 

r**** ' in a fixed extra dimensional background. This extra dimensional background is commonly assumed to be either a 

Minkowski space-time or an anti de Sitter Universe (for another type of background see [9( for instance) . These two 

/~~ ■ cases receive strong motivation from either supergravity or string theory in more than four dimensions. When the 

i | brane is not taken as a probe brane, the extra dimensional geometry can be influenced by the matter fields on the 

CO . brane and a back-reaction must be taken into account. In five dimensions, these effects are taken into account by 

imposing the appropriate Israel boundary condition |6j. On the other hand, for branes with a low enough tension, 

and/or no matter on their world surface, the branes whose dynamics are governed by their Dirac-Born-Infeld action 

evolve in the extra dimensional background leading to a mirage cosmology on the brane [lfj, [ll|. In particular, the 

cosmology on the brane defined by the induced metric can be non-trivial, e.g. inflation can result for instance. 

Recently, the discovery of the acceleration of the Universe [12[ has led to a reappraisal of the underlying assumptions 
of cosmology. In particular, the cosmological principle, which states that the Universe is homogeneous and isotropic 
on large scales, has been under scrutiny [l3|. As a result, Tolman-Bondi models of the Universe where the earth would 
lie close to the centre of an inhomogeneous Universe have been considered jLll - llq ]. Along these lines, it seems to 
be timely to question the usual hypothesis that extra dimensions should be a symmetric space. In 5d, homogeneous 
and anisotropic brane models have been extensively studied [17h20| . in particular with the aim of extending Wald's 
theorem which states that a space-time with a positive cosmological constant and a matter density satisfying the 
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strong energy condition always isotropises [21|] . Other approaches such as in [22[ tackle the issue of isotropisation of 
three large dimensions in the brane gas cosmology framework by allowing an arbitrary amount of initial anisotropy. 
There it was found that the anisotropy reaches a maximum early in the evolution and then approaches zero at later 
times. 

Here and in the following, we will consider that the Universe is a 3-brane embedded in an extra dimensional space 
with no homogeneity at all. We study the dynamics of the brane as it responds to the extra-dimensional cosmology. 
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To simplify matters, we assume that the (d + 4) dimensional space-time is a matter dominated, inhomogeneous and 
isotropic Tolman Bondi Universe with a special role played by the origin of space. On the brane, the induced metric is 
both inhomogeneous and anisotropic. In principle, we should also introduce matter fields on the brane and study the 
back reaction problem. As a first step, we will assume that the brane is a probe and as a result we only investigate the 
mirage cosmology on the brane. It turns out that the dynamics are governed by non-linear partial differential equations 
whose long time behaviour can be analysed in the vicinity of the origin. We find that the asymptotic behaviour of the 
induced metric around the origin becomes isotropic and homogeneous. This isotropisation on the brane is in fact an 
attractor behaviour. 

Of course our model does not directly apply to our 4d Universe. Strong constraints on isotropy during Big Bang 
Nucleosynthesis should be taken into account [2J, |2J|. Moreover, the dynamics would be certainly changed by the 
presence of matter on the brane. Nevertheless, the fact that an isotropic and homogeneous space-time emerging from 
an initially inhomogeneous and anisotropic Universe around the origin of space is certainly appealing. 

The outline of the paper is as follows. In a first part, we describe branes in a Tolman Bondi space-time. In section 
III, we study the mirage cosmology of a neutral brane in a Tolman Bondi Universe. In section IV, we simplify the 
dynamics and focus on the vicinity of the origin of space. In this case, the resulting inhomogeneous and anisotropic 
cosmology converges to a 4d Fricdmann-Lemaitre Universe. Our results are confirmed numerically 

II. BRANES IN TOLMAN-BONDI 
A. Brane Dynamics 

We consider the dynamics of a single uncharged 3-brane described by its world- volume action 

S b = -T 3 j ' d A x^g~, (1) 

where T% is the brane tension and g^ v is the induced metric on the brane related to the extra dimensional metric Gab 
by 

dX A dX B 
9^-G AB — — (2) 

where x M are the brane coordinates and X (a; M ) the brane embedding. Previous studies have mainly concentrated 
on the embedding of 3-branes in Minkowski or anti de Sitter space-times. Here we concentrate on the possibility 
that matter may exist in the extra-dimensions and may lead to a dynamical time evolution of the extra-dimensional 
geometry In the extra-dimensions (the bulk), the action is simply 



Si+d — 7r~2 — / d Xy—GRd+i + S™+4 cr . (3) 
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d+4 

Notice the dimension of K d+4 is — (d + 2). The matter action may comprise ordinary matter like a pressure- less fluid. 
We will also assume that the distribution of matter in the bulk could be inhomogeneous, i.e. the cosmological principle 
could be violated in the bulk. If this is the case, then the brane geometry will be both anisotropic and inhomogeneous. 
Even if the cosmological principle is satisfied in the bulk, we will find that the brane geometry is gencrically anisotropic. 
In order to retrieve gravity in the 4d sense at least in the low energy-long distance regime in a way akin to the 
Dvali-Gabadadze-Porrati (DGP) model |25|], one may include an Einstein- Hilbert term on the brane and couple 4d 
gravity to matter on the brane 

1 



d^xyZ-gRi + s? atter , (4) 

where R4 is the Ricci scalar of the induced metric on the brane and K4 the 4d gravitational constant. Here g™ atter is 
the contribution from the matter species living on the brane. 

In the following we will treat the simpler case where the effects of brane gravity and matter are neglected. In this 
mirage cosmology setting where the dynamics on the 3-brane are simply due to its probe nature and its motion in 
the bulk, we will be able to analyse interesting phenomena such as the isotropisation of the 4-d geometry. Of course, 
the fact that we do not take into account S4 means that our model cannot describe the physics of the Universe since 
Big Bang Nucleosynthesis. To do so, a full analysis of the model involving the total action 

St = Sd+4 + Sb + Si, (5) 



must be tackled. This interesting possibility is left for future work. In the following, we only consider the dynamics 
driven by the brane action Si,. 



B. Tolman-Bondi Space-Times 



Tolman-Bondi space-times in extra-dimensions have been extensively studied |26H32| . Here we review some of their 
salient properties. We consider a (<i + 4)-dimcnsional universe filled with a pressurclcss fluid and a spatially dependent 
energy density p(r, i) whose metric is given by the extension of the Tolman-Bondi solution to (d + 4)-dimcnsions 



ds l 



-dt 2 + e x(r > l) dr 2 + R 2 (r, t)dfL 2 d+2 , 



(6) 



where dVi 2 d+2 is the metric on the (d + 2)-sphere, i.e there are n = d + 2 angular coordinates. The Einstein equations 
yield 



e x = 



R' 2 



l + /(r)' 

where ' = dj dr. The time-dependence of the metric is specified by the Fricdmann equation 



#-m+m 



R 



71—1 



(7) 



(8) 



where ' = d/dt and where f(r) is an arbitrary dimensionlcss function of the radius r. When f(r) = 0, the Tolman- 
Bondi Universe reduces to a Friedmann-Lcmaitre model with no inhomogencity if we assume a cosmically simultaneous 
Big-Bang time, i.e. t = t (r) ^constant independently of r. The Big-Bang time appears as an integration function 
to = to(r) so that each shell has its own big-bang time, a feature which differs from what happens in the Friedmann- 
Lemaitre limit unless to is explicitly chosen to be r- independent. In this case we have R(r,t) = ra(t) where a(t) is the 
scale factor of the Friedmann-Lemaitre Universe. 
The mass function M(r) is defined by 

nM' 



2Kl +d RPR" 



or equivalcntly 



M (r) = 2 -^±± / 



(9) 



(10) 



where p(r,t) has dimension (d + 4). In the flat Friedmann-Lemaitre case, we have p(r,t) — p(t) oc a ("+ 1 ' and 
M(r) oc r n+1 . We can integrate the Friedmann equation to obtain 



„/ n / n + 1 M , 



2/(n+l) 



(11) 



such that at t = to, R(to) = r. Choosing the time origin to be t = 2y t /r n+1 /M/(n + 1), we can write 



R(t,r) 



to 



2/(n+l) 



n + 1 



2/(n+l) 



M l/(«+l) f 2/(«+l)_ 



We then obtain 



,iy /(n+i) 



(12) 



(13) 



which specifies the homogeneous Friedmann-Lemaitre models. 

When f(r) =^ 0, we can follow the approach of Dcbnath and Chakraborty [33( where the inverse of the Fricdmann 
equation |(HJ) is used to obtain t = t(R, r) 



t-tn 



= 2 M-^ r (n+l)/2 ^ 

n + 1 


"1 1 + n 3n - 1 

2' 2(n- 1)' 2(n-l)' 


M 


2 M -l/2^(n+l)/2 p 

n+1 


1 1 + n 3n - 1 
2'2(n-l)'2(n-l)' 


R n ~ 1 f 
M 



(14) 



where at t = to, R(to) = r. We can now simply use R' — (dt/dr)/(dt/dR) to obtain the spatial derivative of R 

3n-l R n - l f 






1. 



n-l'2(n-l) : 



M 



M V R n -1 



^ f M' 

1+71 — -2 — 
^ ' f M 



from which we can obtain A. When n = 3 and therefore d = 1, (|14[) simplifies considerably and we obtain 



R 2 = (t-t Yf±2(t-t )VM 



(15) 



(16) 



which gives an explicit description of the five-dimensional geometry. When / = 0, we find that R 2 (r,t) = 2(i — t a )\/M 
corresponding to a scale factor growing in (t — t^) 1 / 2 . 

In the following, we will consider models where f(r) vanishes at the origin of space r = in such a way that 
f(r) <C M (r) / R n ~ 1 close to the origin. As M{r) oc r n+1 there, this implies that f(r) must go to zero faster than r 2 
locally. This guarantees that the curvature vanishes locally in a neighbourhood of the origin and the extra-dimensional 
geometry is locally the one of a flat Friedmann-Lemaitre model. In particular, this choice specifies that the origin, 
on top of being the centre of the distribution of matter, is a special point of space where the geometry is nearly 
homogeneous and isotropic. The geometry at large far away from the origin remains isotropic but deviates from 
homogeneity as soon as f(r) becomes relevant in the Fricdmann equation. 

III. MIRAGE COSMOLOGY 

In this section we will construct the mirage cosmology set up, first describing the higher dimensional embedding of 
our 3-brane and second, by deriving the dynamical equations that rule the induced 4 dimensional metric. 



A. Brane embedding 



We consider a D3 brane embedded in the Tolman-Bondi background in (d + 4) dimensions that we have presented 
in the last section. The extra dimensional matter distribution is comoving with the energy-momentum tensor 



A„.B 



T A » = p(r,t)u A u», 



(17) 



where matter particles are comoving in the bulk and spherically distributed with u — Sq. In order to preserve 
spherical symmetry of the induced matter distribution on the brane with the induced energy momentum tensor 

dX A dX B 



T^ = Tab 



dx^ dx v 



(18) 



where X A {x^) is the embedding of the brane with coordinates x^ in the Tolman-Bondi Universe with coordinates 
X A , we choose to embed the brane as 



x° = t, 



where x° and x 1 are the brane coordinates. Matter is then spherically symmetric on the brane with 

T^ = p(x°,x 1 )u ll u v , 



where 



u A 



8X A 
dx^ 



(19) 
(20) 

(21) 
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is explicitly u^ = —5® as ua = —S A and dX°/dx ,J ' 

This choice of embedding implies that the 3-brane corresponds to a two dimensional cross-section of the (d + 2)- 
sphcre. Using the (d + 2)-rotational invariance, we can always rotate the cross-section to align it with a planar section 
of the (d + 2) sphere paramcterised by the angles 8 1 and 9 2 . Figure Q] illustrates this embedding. Our 3-brane consists 
of the cone defined by the angles 8 3 , ..., 9 n = constant and the free angles 8 1 and 9 2 . 




FIG. 1: Schematic description of our higher dimensional embedding. Our 3-brane is represented by the cone defined by the 
fixed angles # 3 , ..., 6 n and the freely rotating 8 1 and 6 2 . 



Denoting by x an d </> the angles of the 2-sphere on the world- volume of the brane at constant x , we find that the 
embedding is characterised by the two fields 



e 1 = 6\t,r,4>,x), 9 2 ^9 2 (t,r, ( p,x), 

The other coordinates on the (d + 2)-sphere are in the normal space to the brane and such that 

e i = eL i> 3, 



(22) 



(23) 



where 0$ are constant angles. 

We are interested in the geometry seen by a comoving observer on the brane characterised by the velocity vector 
i/' = <5g which differs from the velocity vector of matter particles as seen from the brane m m = <?'*<5q when 6* 1,2 are 
not time independent. The spatial geometry seen by such a comoving observer is defined by the metric 



hf = g" v + v»v v , 



(24) 



satisfying h^ 



aauioiying n^v = 0. On the space orthogonal to v^ it reduces to g lJ where x 1 ' 3 = (r, x, 4 1 )- This metric is both 
inhomogeneous and anisotropic. Moreover proper time for a comoving observer at rest on the brane 



dr l 



-gudt 



(25) 



is also defined inhomogeneously in space. This is due to the (t,r) dependence of 9 1 ' 2 . 

In the following, we will concentrate on the vicinity of the origin of space where we have assumed that the curvature 
effects in the bulk can be neglected. In the neighbourhood of this point, the bulk Tolman-Bondi metric reduces to a 
Fricdmann-Lemaitre space-time. On the other hand, on the brane, the local geometry defined by the induced metric 
is neither homogeneous nor isotropic there. An observer analysing the geometry of its local patch would see the 
Universe as both inhomogeneous and anisotropic. This is the case despite the fact that matter as seen on the brane 
is spherically symmetrically distributed according to p(r,t). We will study the long term dynamics of the geometry 
around the origin and see when the local geometry isotropises and becomes homogeneous up to points where the 
curvature effects in the bulk measured by f(r) cannot be neglected anymore. 

The dynamics of the brane, solely governed by its world volume action and depending on the tension T3, depends 
on the induced metric. With the chosen embedding, we have for the induced metric: 



ds = gudt + ijrrdr + 2g tr dtdr + g^d<\> 
9xx d X 2 + ^94>x d $ d Xi 



(26) 



which can be related to the (4 + d) dimensional metric Gab using ([2]) . This gives explicitly 

2 / d6 ? x 2 

dr dr \ dr J \ dr 

where g a b is the metric on the {d + 2)-sphere and the 9 a, s are spherical coordinates. Here we have the diagonal metric 
with gn = 1, gi2 = and 522 = sin 8 1 . Similarly we have 



g rr = e^ + R 2 g ab — — = e x + R 2 { — ) +R2 s in 2 e 1 [ — ) , (27) 



d° a de» _ 1 , d2 /^V , D 2„;2 Q ifde 2 



„» a6» a a6< 6 „ 2 d0 1 d$ 1 „- . 2n ,d9 2 d6 2 
at ar ai ar at ar 



and 



^^^^^f^y+^sin^if^V, (30) 



0*»a^ = fi2 /90iy fi?2sin20l W 

% ax \9x/ \<9x 

, 2 a6» A a6» s .a^afl 1 2 2 , a6» 2 a6> 2 



^ = R 2 ( d -f] + R2 ^ 2el ( d -f) ^ (31) 



With these expressions, we can express the volume element of the brane in terms of the fields 9 12 . The dynamics of 
the brane is simply obtained by looking at extremal world- volumes. 

B. Mirage dynamics 

We are now ready to write down the 3-brane world- volume action (TTJ , find the corresponding Lagrangian in terms 
of the higher dimensional angular coordinates and compute their equations of motion. The action of the brane is given 
by the integral of the volume element 

V^ = [-(9tt9rr - 9tr)(94"p9xx ~ ff$ X )] ) ( 33 ) 

leading to a Lagrangian 

£ = -T3 [-{gttgrr-g 2 r){94>4>9xx-9lx)] • ( 34 ) 

We expand this Lagrangian to second order in the time and radial derivatives as we are interested in low energy 
modes. Hence we get the Lagrangian 

C = -T 3 i?V/>in^||l + ^(^) 2 + ^sin 2 ^(^ " 

1 u2(de 1 \ 2 i 2 . 2o if9e 2 \ 2 \ (de l de 2 de 1 de 2 \ 
-2 R {-at) - 2 R sm e {-sT ) ) [-84 fl* - -^WJ ' (35) 

which can be conveniently cast in the form 

C = -JKL, (36) 

where J captures the angular properties of the embedding 

J -{d^ d x ~ d x d<t>)' (37) 

while K/T3 is the volume element of an isotropic and homogeneous Friedmann-Lcmaitrc Universe in 4 dimensions 

K = T 3 R 2 e x/2 s sin 6 1 , (38) 



where s — | sin 9 X \/ sin 9 1 = ±1, and finally we have a reduced Lagrangian L up to second order in the time and radial 
derivatives which encompasses the dynamical properties of the 3-brane 

The dynamics of the brane are governed by the Euler-Lagrangc equation for 9 1 and 9 2 . Noticing that the Lagrangian 
does not depend on 9 , we find a conservation equation 

d^A" = 0, (40) 

where 

A» = 9C (41) 

is the conjugate vector to 9 2 . We will see that this vector is not automatically a constant, it has an explicit space-time 
dependence gencrically. Moreover, different choices for A* 1 would lead to different solutions to the brane dynamics. In 
the following, we will see that the dynamical equations can be solved for a subclass of vectors A^ . The general case 
is certainly interesting although challenging technically. Our analysis is general in this section. Explicitly, we have 

(42) 

(43) 

(44) 

(45) 



A 1 -- 


= JKR 2 sm 2 9 1 ^—, 
dt 


A r -- 


= -JKR 2 e~ x sin 2 O 1 ^- 


A^ -- 


ox 


A* -- 


= —KL——. 



Equation (f3"?| can now be written more compactly using (|4"4")l and (|4"5j) . 

KL v "< x 



J=--L(A*9 2 X + AH 2 ^). (46) 



and then (|4"2")l and (|4"3"j) can be recast as 

a, LA* 



R 2 sin 2 9 1 {A*9 2 



x 
LA r 



R 2 e- X sin 2 9 1 {A^9 2 X + A^9\) 



(47) 
(48) 



which upon substitution into Eq. (|39[) give a quadratic equation for L whose solutions are 



L = (1 ± VI - 4AC*)/2A, (49) 
where we have conveniently defined 

A = eA (^) 2 -^) 2 m 

2R 2 sin 2 9 1 (Ax9 2 x + A^9 2 (f> ) 2 ' ( ' 



C = l + \R 2 



,-A//il \2 //il\2 



r-(o i r . (si) 



In the following, we will retain the + solution and impose that \AC\ <C 1 is such a way that L ~ 1/A. The Euler- 
Lagrange equation for 9 1 is a very complex partial derivative equation which does not seem to be tractable in full 
generality, however, by performing suitable simplifications we will find an elegant expression capable of being studied. 



C. Simplified dynamics 

We will reduce the dynamics of the brane in a Tolman-Bondi bulk by using simplifying assumptions. The first step 
consists in choosing a subclass of vectors A^ for which the dynamics can be solved in a non-trivial way. So we will set 

A r = A x = (52) 

which immediately implies that 

8 2 r = B\ = 0, (53) 

respectively and 

A= -(A t ) 2 /2R 2 sm 2 6 1 {A< t ') 2 (e 2 <j> ) 2 . (54) 

In addition, we will assume that e~ x l 2 Q\. <C 9 1 . In effect, we consider that the dynamics is very close to being 
homogeneous. These choices will allow us to simplify the equation of motion for 9 1 , which can now take the form 



i R + 2 x+2 w + ¥ir + ^^ {9) ~ e sm e ¥-m [w: • (55) 



Notice that the radial partial derivatives are not present. When \AC\ <lwe find that 9 1 can be approximated by 



3l 



(A*) 2 



• x 2T 3 A$eV 2 R 4 s sin 3 ^ 1 ^) 2 ' 

In the following we will assume that A* is constant in time. Differentiating 9 with respect to time we find that 

Q2 



(56) 



*-- * [ 4 r + 2 +3 ^¥ + 2 ¥1 + M ' (57) 



which upon substitution in ([55]) gives 

9 -{ 4 R + 2 X + 3 9% +2 A^) e ~ 5 ^¥ {6) ~ 6 Sm ° e^diyWj- 

This is the partial differential equation which governs the brane dynamics. 

We have seen in Eq. (|47l) that 8 2 = —LA t /A (f '9 2 ,R 2 sin 2 6 1 , which under our approximations gives 



(58) 






~ 2 ~it e %- (59) 



As we have required 8\ = we find using Eq. ([56]) that 9 2 must be linear in <j), i.e. 6 2 i = cq where cq is a constant 
and 

9 2 = cocj, + 2c / -rrdt. (60) 

J to A 

Notice that Q\ = implies that A^ j A 1 must be r-independent. 

We have now obtained the equations of motion for 9 1 = 1 (t, x) an d the general form of the solutions for 2 = 
2 (t,<fi,x)- One must therefore solve 



'£ + i A+a £r -«=>>■- 



2/jl J °,X 



,X 



The source term depends on A^ and A t . We have therefore reduced the brane dynamics to a single differential equation. 



IV. THE BRANE GEOMETRY CLOSE TO THE ORIGIN 

In this section we will solve (|61|) by decomposing 9 1 and A^ in their respective time and x-dependent parts. We 
will see that the solution has a late time attractor. 

A. Local Attractor 

Even the simplified dynamics are difficult to analyse. We concentrate on the interesting case of a 3-brane embedding 
in a Tolman-Bondi Universe where the curvature effects measured by f(r) can be neglected in the neighbourhood of 
the origin. In this case, the extra-dimensional dynamics in a local patch becomes: 

t \2/(«+l) , 4 x4/(n+l) 



*<*' r > = r UJ ' e UJ ' (62) 

and we focus on the induced metric on the brane which takes the simplified form 

gtt = -1 + R 2 (0 1 ) 2 + R 2 sin 2 e 1 ^ 2 ) 2 , (63) 

9rr = e\ (64) 

g H = R 2 sm 2 d 1 (6 2 4> ) 2 , (65) 

~9xx = R 2 {e] x ) 2 + R 2 *in 2 6 l (el)\ (66) 

9lPx = R 2 sm 2 e 1 e 2 x e%. (67) 

We will analyse the geometry on the brane when we have the explicit solutions of the equations of motion. 

We are interested in embeddings where the lapse function does not deviate too much from 1. We consider then 

where the forcing term is 



-,2 nl iX 



We look for a solution for 8 1 locally such that 9 1 ta do and we write 



2T 3 s sin 3 eoM^cl^ge^R 4 J f 



(69) 



d X + ait), (70) 



where #o is a fixed angle arising as an integrating constant, a(t) an unknown function yet to be determined and we 
have introduced 

A*(r,t, X )=M<>r i f(x)g(t), (71) 

in order to separate the time and x - dependent contributions that make up A?. Here, Mo is an unspecified energy 
scale. Notice that the r-dependence cancels if and only if A f varies in r 4 as required by 6 2 r =0. We now consider that 
X measures the small deviations from #o and approximate / = 1 + fix, 

gl=fl0 - 2T,^ a ^ r4geV2jz4 X + a(t). (72) 

We identify the constant of motion A* as 

(A 1 ) 2 = 2T 3 Aiy o s sin 3 e r 8 g , (73) 

and write the solution for 6 1 as 

e'=9 - 9 je^/ 2 (^y x + a(t), (74) 
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which is r-indcpcndcnt, hence satisfying our assumption on the small r variations of 9 1 . At t — <o wc have that 



1 = - X, (75) 

2 = c o 0, (76) 

where we have chosen a(to) = 0. Thus initially, the spatial part of the induced metric is 

9rr = 1, (77) 

~g H = ^sin 2 ^, (78) 

9xx = r\ (79) 

.90x = 0, (80) 

which corresponds to flat space. As time passes, the metric is deformed and deviates from fiat space. On the other 
hand the lapse function g t t can not be uniformly equal to —1 even initially, implying that time flows differently at 
different points of space. 

The analysis of the equation of motion for 6 l is greatly simplified if one studies the evolution of small deviations 
from #o = 7r/2, i.e. we consider an expansion in small \ around 9 1 = ir/2. The general case around a general 6q follows 
easily as we will see soon. Substituting the previous ansatz for 9 1 in the equation of motion and separating the terms 
in x from the ones of order zero in x, we obtain two second order differential equations for a and g 

a- U| + iA + 2^)d = 7 Q , (81) 

where I a = -(2f 1 g Q M^/T 3 )e^ 2 {g/g ) 3 (R/r^ and 

g\ 2 k JrY lv fq A R 1;V 

- +4 4 — +-A--+4— + -A - 

9 J R \R 2 \g R 2 j 



h A R I- \ g R 1; 

- + 4 h -A - + 4 h -> 

9 R 2 ) \g R 2 



2 -+ 4 15+oM (:+4- + -A|=/ 9 , (82) 



where I g = —(4:figoMQ/T 3 )e x (g/go) 4 (R/r) s . The solutions to these equations can be easily obtained. 

First of all it follows from the a-equation that the homogeneous solution when the forcing term vanishes is 

a{t) = do f dt'e x ' 2 (-} (lL\ . (83) 

The ^-equation (1821) is equivalent to a Riccati equation 

V = 3irV A/ V 2 + R 4 e x/2 I gi (84) 

and we have defined 

y = i?4eA/2 (f +4 | + ^)' (85) 

Putting b = 3i? _4 e _A / 2 and c = R 4 e x ^ 2 I g , the Riccati equation becomes a Sturm-Liouville equation 



X 

where the following correspondence 



be -*/'*-( J- 



X = 0, (86) 



u,' 



V = ~, X = b-^oj, (87) 

bw 



applies. The homogeneous solution (c = 0) is 



X = b- 1 ' 2 ( ci + c 2 / dt'b{t') , (88) 
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where c\ and ci are integration constants. As a result, we find that the solution of the homogeneous g-cquation is 

-1/3 



90 



. x 10/(n+l) 
l \ 



1 + C 3 



. N (9-n)/(n+l) 



1 



(89) 



where c^ is a constant. We find that as long as the term within the power —1/3 does not vanish, a statement which 
only depends on the initial conditions, the asymptotic behaviour of the homogeneous solution is 

g~ lg t- w l^\ (90) 

where 7 9 is a constant. Using this result and substituting in (|83j) one finds that the asymptotic behaviour of a is 

(91) 



a « 7a <("- 9 )/(" +1 > + , 



where 7 Q and a^ are constants, i.e., for n < 9, a converges to a constant at infinity. As I g = 0(t 20 /( n + 1 )) falls off 
faster than i~ 2 and I a = O(t~ 20 ^ n+1 ^) faster than a = 0(H™ +11 )/(" +1 )) for n < 9, we find that the homogeneous 
solutions (|90|) and (|9T|) are late time attractors where the influence of the forcing term I g and I a becomes negligible. 
The same results hold for any 8q as wc can neglect the quadratic term in Q in Eq. (|6Tj) as the solution of the 
homogeneous equation behaves as 6 1 = 0{f5\ + / 5 2 t(' R_9 )/(' re + 1 )) while / = O(t~ 20 /( n+1 ' ) ). Hence, in the neighbourhood 
of each point 9q, we find that 8 1 and 9 2 become time-independent. As a result the metric converges to 



dsi 



-dt 2 + e x dr 2 + R 2 [sin 2 ^cgtty 2 + ((0 1 ) 



and on the attractor B 2 evolves as 



X't^-oo 



>*(f 



1/2 



,2/)l/ 



)Loo) d X 2 + co sin 2 Wjt-.ood^dx] , (92) 



.X'i-^oo 



(9-n)/(n+l) 



(93) 



which indeed tends asymptotically to a constant value (9 2 )t_>oo- Also 6*1. is equal to minus one on the attractor. As 



a result the mapping between (x, </>) and 



') becomes in the late time regime 



9 1 

9 2 



X + Q'o 



co^+(0f x )t->ooX- (94) 

This is a linear transformation with constant coefficients. The transformation of the metric g a b on the two sphere 



ft3-5o6 fel Qxj 

where x 1 = (0, x) is a linear transformation defined globally for any time. The asymptotic metric is then 

ds^ = -dt 2 + e A dr 2 + R 2 g ij dx i dx :i , 
which is equivalent to 

ds 2 ^ = -dt 2 + e x dr 2 + R 2 g ab d0UoodOUoo- 



(95) 



(96) 



(97) 



Hence we find that the asymptotic metric is isotropic, the induced space-time on the brane isotropises and converges 
to a Friedmann-Lemaitre space-time close to the origin, of the radiation dominated type for n = 3. 

We will see in the next section that the solution obtained here locally around each point specified by 6q in a small 
X expansion is in fact a global solution of the equations of motion. The analysis of this section shows that the global 
solution wc will find in the next subsection is an attractor. 



Global Attractor 



We have seen that the fields 8 1 and 8 2 become locally time-independent in any neighbourhood of #o and that this 
solution is an attractor, i.e. all solutions around #o converge towards a unique long time solution determined essentially 
by the asymptotic behaviour of the function g(t) ~ i _10 /(™+ 1 ). Here we will find a global solution which corresponds 
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to this attractor, i.e. we will find a solution for very large time which matches the local behaviour of the solution 
determined by g(t) ~ £- 10 /(™+ 1 ). I n particular, we will find that the induced metric on this attractor is the one of a 
radiation dominated Fricdmann-Lcmaitre Universe when n = 3. 

The homogeneous solution of (|5T|) when the term in (9 ) 2 is neglected is simply 

d\t, r, x) = E( X ) f dt'R^'^A^f + F( X ). (98) 

The first term is time dependent while the second one is only a function of x- This general solution must be compatible 
with the differential equation ([55)1 which can be obtained putting E(x) = 0. In this case, the behaviour of 6 1 (x), 
necessitates that (A*) 2 jA^ e x l 2 R 4 (9 2 ,) 2 be time and r independent. This is achieved by choosing 9 2 < = cq constant 
and A 1 time independent. Moreover we choose A 1 ex r 4 and A^ oc r 4 implying that 9\^ is r- independent. Similarly A^ 
must be such that A^e x ^ 2 R 4 is time-independent. This implies that 

^ocr 10/(n+1) . (99) 

Notice that the behaviours of A^ and A* as a function of r and t match the asymptotic behaviour of the local solution 
as they are both proportional to r A and that A^ has the same time dependence as the asymptotic behaviour of g(t) 
around any background value 9q. We can now integrate exactly ((55)) and get the global mapping between x an d 9 1 

(A*) 2 
d9 l 1 sin 9 1 1 3 = T-rr-rrr-rX, (100) 

where the right hand side is independent of r and t. Expanding the left hand side around a given 9q to linear order 
gives the local solution of the previous section. 

For this solution, the non-linear term in (9 ) 2 of (|55|) vanishes identically while the forcing term converges to zero 
/ = 0(t -2O /( n+1 )). Hence this global solution is valid at infinity and around each 9q it is identical to the one obtained 
in the previous subsection from (|89[) with C3 = 0. Hence, we have found a global solution of the equation of motion 
valid in the late time regime and we have shown that this solution is an attractor as the solution (|89p converges to 
(OH). 

In conclusion, we have found that locally around the origin of space-time on the brane, the geometry isotropises in 
the long time regime. This implies that a local observer, comoving on the brane would eventually see a flat Friedmann- 
Lemaitre space-time in its local patch. This result is special to the origin where we have assumed that the bulk effects 
of the curvature of space-time can be neglected. In a sense, we have shown that for a brane on which the matter 
density induced from the bulk is spherically symmetric and where the effects of the bulk curvature are negligible, the 
brane geometry eventually becomes homogeneous and isotropic locally. Of course, this is not the case away from the 
origin implying that interesting inhomogeneous effects on the propagation of light on the brane for instance may be 
induced. The study of the brane geometry at large away from the origin is left for future work. 

C. Numerical Results 

In this section we will illustrate the previous results and show numerically that the attractor mechanism is at play 
for a wide range of initially anisotropic and inhomogeneous space-times. First, wc must demand that the late time 
solution is such that g tt « -1, i.e., that R 2 {9 1 ) 2 < 1 and R 2 (9 2 ) 2 < 1. As we have seen, A 4, = M e r 4 g(t)f(x) and 
(A*) 2 ra 2T 3 M 6 c 2 l r 8 # then we have from ||S5J) that 

Of2 71 f 6 / 4. \ 16/(n+l) 

W-^ftr'lT «!■ (101) 



T 3 au V * , 

Second, we must also require 4|AC| <C 1 so that we can use the approximation L sa \j A and therefore, we have the 
additional constraint 

4T 3C 2 1 / t \W(»+l) 



41^1 « 77^— H «!• (102) 

M°f 2 r 2 g \toJ 

When numerically solving the equations of motion we have started the evolution at £0 = W^^t^BN where the time 
of Big Bang Nuclcosythesis is tgsN ~ 10 26 GcV~ 1 . We will compute solutions for the scale of the Hubble horizon at 
that time, i.e. r = to. We set T 3 = lGcV, Mq = lGcV, n = 3, /1 = 10~ 15 , x — and in order to satisfy conditions 
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,x 10' 




10"' 

t/tBBN 



10" 



10" 



FIG. 2: Evolution of a{t) with time for T 3 = lGeV, M = lGeV, n = 3, /i = 
ff(t = 0) = -9.5g /(n+l), c = 4.4 x 10" 8 , a(t = 0) =0, a(t = 0) = 10~ 4 and r 
a constant value at late times as we expected from Eq. (|9ip . 



10"^, X = 0, go = g(t = 0) = 1.2 x 10" 4 \ 
= 10 4 tBBN- One observes that a approaches 




FIG. 3: Evolution of g(t) for the same parameters of Fig. [2] We confirm that g(t) decays as a power law given by Eq. (f90 



(jlOip and (|102|) we have used go = 1-2 x 10 -41 and Co = 4.4 x 10 -8 . The small size of /i is chosen to prevent the 

source terms I a and I g to be of any importance. As initial conditions we have used a(t = 0) = 0, a(t = 0) = 10~ 4 , 

g(t = 0) = go. If the evolution was already on the attractor we would have g = —10go/(n + 1). Instead we start 

slightly away from the attractor solution using g(t = 0) = — 9.5go/(n + 1). The result of the numerical solutions is 

shown in Figs. [2] to [5] In Fig. [5] we can observe the time dependence of a(t) and confirm that it quickly approaches a 

constant value a^ as we expect from Eq. (pJTj) . On the other hand, it can be seen in Fig. [3] that g(t) approaches the 

power law solution g(t) ~ t~ 5 ^ 2 as we have seen in (|90|) for n = 3. 

Figure |4] shows that 1 never deviates substantially from 9q = 7r/2 which ensures that our approximations, i.e. 

neglecting the (9 1 ) 2 term in Eq. (|6Tj) and making sin 3 9 1 ps sin 3 #0, are valid and under control. 

We also demanded that 4|AC| <C 1 so that L rj 1/A which allows us to have a simple expression for 9 1 and be able 

to write Eq. (|56p. As we can see in Fig. [5] this is satisfied for our choice of parameters. 

Figures [51 [7] and [5] show that indeed the metric on our 3-branc quickly isotropises at late times as gtt —t — 1, 9xx/9<P4> ~^ 

constant, g$ x lg§§ — > constant and g^, x jg^ too. This is what allows us to redefine the angular variables asymptotically 

and find that the spatial metric is the one of flat space. 
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t/tBBN 



10" 



10" 



FIG. 4: Evolution of A = 100(S 1 — n/2)/(n/2) for \ = which is the deviation of 9 from 9q = 7r/2 in percentage. This example 
clearly shows that the deviation is very small which ensures that our approximation sinfl 1 » sin So is valid during the whole 
time interval analysed. 




FIG. 5: The time evolution of 4|AC| which we require to be smaller than unity for our simplification L ~ 1/A to be valid. This 
is indeed the case for the time interval we are covering and for our choice of parameter values. 



CONCLUSION 



We have studied the dynamics of a 3-brane embedded in an inhomogeneous extra dimensional Universe. The mirage 
cosmology on the brane reveals that the late time behaviour metric on the brane becomes isotropic and homogeneous 
in the vicinity of the origin, where the effects of curvature in the bulk can be neglected, even though the initial 
condition were not so. This isotropisation process is extremely interesting in view of the possible existence of extra 
dimension and the strongly motivated possibility that matter, even cosmological matter, could exist non-only in 4d 
but also in the extra dimensions. In this case, the dynamical interplay between the brane degrees of freedom and 
the extra dimensions becomes more complex than in the usual case of an extra dimensional symmetric space, e.g. 
anti de Sitter or Minkowski. Indeed, one may envisage situations where both matter would exist on the brane and 
in the extra dimensions, with possible exchanges between them. In this paper we have only analysed the mirage case 
where matter on the brane is subdominant and the dynamics are governed by the embedding of the brane in the 
time-dependent extra dimensional background. The coupling between matter on the brane and outside the brane in 
an inhomogeneous context is left for future work. The presence of inhomogencities in the extra dimensions may even 
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FIG. 6: Evolution of the metric component gtt- To recover a Friedmann-Lemaitre metric this quantity must be —1 which indeed 
happens for our solution at late times. 




10"* 10 10 10 10 

t/tBBN 

FIG. 7: The time evolution of g X x/9<t><t>- This ratio approaches a constant at late times as expected from Eq. (|92 



shed some light on thorny issues such as the acceleration of expansion of the 4d Universe or the role of dark matter. 
These tantalising possibilities are beyond the present work. 
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